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DETERMINATION  OF  OPTICAL  CONSTANTS  FROM  EXTINCTION  MEASUREMENTS 


1.  INTRODUCTION 

The  Interaction  of  electromagnetic  radiation  with  spherical  particles  is 
described  by  the  well-known  Lorenz-Mi®  formalism. The  optical  properties  of  the 
particulate  material  enter  into  this  formalism  by  way  of  the  optical  constants  (n, 
k)  which  are,  respectively,  the  real  and  imaginary  components  of  the  spectral 
complex  refractive  index. 


N(X)  -  n(X)  -  ik(X) 


(1) 


where  X  denotes  the  wavelength.  With  the  advent  of  modern  high-speed  computers  and 
the  development  of  reliable  Mie  scattering  codes, it  has  become  a  commonplace 
numerical  procedure  to  predict  the  extinction  produced  by  an  ensemble  of  spherical 
aerosol  particles  once  the  optical  constants  (n,k)  and  the  particle  size 
distribution  are  known.  The  inverse  problem  of  determining  the  optical  constants 
from  the  Lorenz-Mle  formalism  when  measured  values  of  the  spectral  extinction  and 
the  particle  size  distribution  are  known  has  proved  to  be  quite  intractable;  this  is 
the  problem  which  this  paper  addresses. 


12-14 

As  noted  by  several  workers,  the  experimental  difficulties 

associated  with  determining  the  optical  constants  by  transmission,  absorption,  and 
reflectance  measurements  performed  on  small  samples  of  the  aerosol  material  lead  to 
highly  questionable  results,  A  few  techniques  which  avoid  these  difficulties  by 
using  the  Lorenz-Mie  equations  to  invert  measurements  of  particulate  optical 
properties  to  find  the  optical  constants  have  been  reported.  Wyatt15  and  Pluchino 
et  al«*6  have  measured  the  angularly  dependent  intensity  of  laser  light  scattered 
from  single  particles  which  were  suspended  in  a  Millikan  oil  drop  apparatus;  the 
Lorenz-Mie  equations  were  then  used  to  fit  the  angular  scattering  data  by 
curve-fitting  techniques  in  which  (n,k)  appear  as  parameters. 


Another  approach  is  to  use  extinction  measurements  which  are 
experimentally  simpler  than  angular  scattering  measurements,  but  they  have  the 
disadvantage  that  the  Lorenz-Mie  inversions  are  not  unique  with  respect  to  (n,k)  as 
they  are  with  the  angular  scattering  measurements,  Janzen1^  has  described  a 

technique  for  determining  the  optical  constants  of  a  carbon  black-water  colloid  by 
using  the  Mie  equations  to  obtain  a  least  squares  fit  to  measured  extinction 
spectra.  Janzen's  technique  overcomes  many  of  the  difficulties  of  the  traditional 
methods,  but  there  are  three  limitations  to  his  approach:  (1)  the  polydispersity  of 
the  particulate  material  is  not  accounted  for,  (2)  the  optical  dispersion  of  the 
material  is  not  taken  into  account,  and  (3)  the  optical  constants  which  are  found  by 
this  technique  satisfy  the  Mie  equations,  but  they  are  not  necessarily  the  unique 
(n,k)  pair  which  characterizes  the  optical  behavior  of  the  material.  In  this  paper, 
we  will  describe  a  method  of  deriving  the  optical  constants  of  an  aerosolized 
material  from  extinction  measurements;  this  method  overcomes  the  limitations  cited 
above. 
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The  extinction  coefficient  (in  / go)  of  an  aerosol  at  wavelength,  X,  Is 


defined  as 


a 


(2) 


where 

2 

s*  •  optical  cross  section  of  the  1th  particle  (m  ) 

•  the  mass  of  the  1th  particle  (gm) 

and  the  sum  extends  over  the  ensenfcle  of  particles  contained  In  the  optical  path. 
For  a  continuous  distribution  of  particles,  the  expression  for  the  extinction 
coefficient  becomes 


a  ■  /  a(z)dM 


(3) 


where 


a(z)  -  extinction  coefficient  for  a  particle  of  size  z 
dM  -  the  mass  distribution  function  of  the  particles 

In  order  to  apply  the  Lorenz-Mle  formalism,  each  aerosol  particle  will  be  considered 
to  be  spherical  with  an  extinction  coefficient  which  Is  given  by 


a 

o 


3Q  [N(X),x] 

_ £ _ 

2pD 


,U) 


where 

Qe  -  the  efficiency  factor  for  extinction 
X  -  xD/X,  the  size  parameter 
p  -  the  particle  density  (gm/cm^) 

D  ■  the  particle  diameter  (urn) 

It  Is  also  convenient  to  define  the  dimensionless  extinction  coefficient 


a.  -  a  e  X  e  p 
D 


(5) 
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The  difficulties  of  using  extinction  measurements  to  determine  the 
optical  constants  are  Illustrated  in  figure  A-l.*  For  a  given  wavelength,  a  value 
of  Op  may  correspond  to  any  (n,k)  pair  lying  on  the  pertinent  isopleth; 
similarly,  for  a  given  particle  else  distribution,  the  (n,k)  pair  corresponding  to  a 
particular  value  of  is  not  unique.  However,  over  a  wavelength  spectrum,  the 
values  of  [n(X),  k(X)]  are  subject  to  the  constraint  imposed  by  the  Kramers-Kronig 
(KK)  dispersion  relation.*  ~  In  the  development  which  follows,  an  iterative 
procedure  will  be  described  which,  after  an  initial  estimate  of  n(X)  over  the 
wavelength  spectrum,  employs  a  simple  numerical  algorithm  to  fit  the  k(X)  spectrum 
to  the  measured  extinction  spectrum  and  the  KK  dispersion  to  establish  a  new 
estimate  of  the  n(X)  spectrum.  It  will  be  shown  that  this  procedure  will  correctly 
determine  the  optical  constants  of  the  airborne  material. 

2.  THEORY 

2.1  Computation  of  n(Xp) . 

The  Kramers-Kronig  relationship  between  n(XQ)  and  the  k(X)  spectrum  is 

given  by: 


n(Xo) 


k(X)dX 

,2 


where  P  Indicates  that  the  Cauchy  principal  value  is  to  be  taken.  The  Integral  in 

equation  6  is  to  be  evaluated  over  the  entire  electromagnetic  spectrum;  and,  since 

k(X)  is  known  only  over  a  finite  region  (X  <  X  <  X  )  of  the  spectrum,  it  is 

min  max 

necessary  to  extend  the  k-spectrum  beyond  the  wavelength  region  where  data  are  known 
experimentally. 

The  inaccuracies  Introduced  by  extending  the  k-spectrum  can  be  reduced  to 
relative  insignificance  by  employing  a  subtractive  Kramers-Kronig  analysis  (SKK). 
Assume  that  the  real  part  of  the  refractive  index,  n(X^),  is  known  at  some 
wavelength,  X^,  and  subtract  the  KK  expression  for  n(Xp  from  the  KK  expression  for 
n(X  )  in  order  to  obtain  the  subtractive  Kramers-Kronig  relationship:^-^1 


»(xo) 


n(Xt) 


Xk(X)dX 


'  2  2' 
Xz  -  X* 

,  o 


The  Integral  is  now  evaluated  by  letting  k(X)  ■  k(X  )  for  0  <  X  <  X  and  k(X)  - 

min  min 

k(X  )  for  X  <  X  <  In  general,  these  are  not  good  physical  approximations 


*  All  figures  are  in  the  appendix. 
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for  k(X)  in  the  spectral  regions  where  k(X)  hes  not  been  measured;  but,  as  mentioned 
previously,  this  approximation  will  have  an  insignificant  effect  on  the  result 

produced  by  the  SR  algorithm, 

2.2  Calculation  of  k(X). 

In  the  previous  section,  it  was  shown  that  SKK  analysis  can  be  used  to 
calculate  n(X  )  once  n(X  )  and  k(X)  are  known.  A  method  of  obtaining  k(X)  from 

O  1 

measurements  of  the  extinction  spectrum,  a(X),  and  the  particle  size  distribution 
will  now  be  developed.  For  a  continuous  distribution  of  spherical  particles,  the 
expression  for  the  extinction  coefficient  of  a  single  spherical  particle  with 
diameter,  D,  (equation  4)  may  be  Integrated  to  find  an  expression  for  the  extinction 
coefficient  of  the  particle  ensemble  at  wavelength  X: 


u(n,k)  -  /  aQ(n,k,D)  dM 


(8) 


For  this  study  the  mass  size  distribution  was  described  by  the  log-normal 
distribution  function 


dM 


2 

-l/2[in(D/Dm)/tno  ] 
e  8 


dinD 


(9) 


where 


Dm  ■  mass  median  diameter  (MMD) 

Og  -  geometric  standard  deviation 

It  will  now  be  assumed  that  an  estimate  of  n(X)  •  n„  is  available  from  the  SKK 
analysis;  and,  therefore,  the  extinction  coefficient  depends  on  k  only 

o(k)  ■  a(nQ,k)  (10) 

The  Isopleths  of  aD,  shown  in  figure  A-l,  depict  the  behavior  of  the  integral 
kernel,  a  ,  in  the  expression  for  a(n,k)  as  the  size  parameter  varies.  Many  of  the 

O  if) 

Interesting  features  of  these  Isopleths  are  due  to  optical  resonance  phenomena. 
The  resonance  due  to  the  first  surface  polarlton  mode  for  X  ♦  0  is  located  at  the 
point  ■  (0,  /F  )  and  moves  according  to  the  equation^ 


k  -  /T(l  +  6/5  X2)1'2  (ID 

to  larger  values  of  k  as  the  size  parameter  increases.  Also,  the  resonances  due  to 
the  bulk  polarlton  modes  can  be  seen  to  move  along  the  n  axis  toward  the  k  axis  as 
the  size  parameter  Increases  (figure  A-l,  E  to  I). 
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Let  us  now  divide  the  n,k  plane  by  a  horizontal  line  through  the  locus 
and  accept  values  of  k  which  lie  In  the  lover  half  of  the  divided  n,k  plane;  i.e.,  0 
<  k  <  /2.  In  this  region,  the  extinction  coefficient  Is  a  single-valued  function  In 
the  vicinity  of  the  first  surface  polariton  resonance  and  no  significant  loss  of 
generality  is  Incurred  by  Imposing  this  condition  since  most  condensed  matter  has  k 
values  which  lie  in  this  region*  If  n(X)  is  now  considered  to  be  fixed  at  some 
value  nQ,  it  can  be  seen  In  figure  A-l  that,  in  general,  the  dimensionless 
extinction  coefficient  Increases  monotonically  with  respect  to  k;  this  fact  forms 
the  basis  of  the  procedure  for  computing  k(X).  Since  the  isopleths  of  represent 
the  Integral  kernel  of  equation  10,  a(k)  itself  is  expected  to  be  monotonic  with 
respect  to  k  over  a  wider  range  of  n  than  is  indicated  in  figure  A-l.  This  has  been 
previously  demonstrated  in  work  by  Jennings  et  al.**  who  computed  isopleths  in  n,k 
space  of  the  extinction  coefficient  at  10.6  pm  for  log-normally  distributed  aerosols 
which  had  number  median  diameters,  D  ,  of  2.0  and  20.0  pm  and  a  «=  1.5.  For  these 
cases,  a(k)  is  shown  to  be  monotonic  for  10" J  <  k  <  / 2,  when  0.1  <  n  <  8  for  Dn  = 
2.0  pm  and  when  0.1  <  n  <  4  for  D  »  20.0  pm. 


The  effects  on  the  monotoniclty  of  the  extinction  coefficient  which 
result  from  changing  the  particle  size  distribution  and  the  real  part  of  the  complex 
refractive  index  are  shown  in  the  curves  of  figure  A-2.  In  these  curves,  the 
dimensionless  extinction  coefficient,  integrated  over  the  indicated  log-normal  size 
distribution  functions,  is  plotted  with  n  fixed  as  a  function  of  the  mass  median 
size  parameter 


X  -  if  D  /X 
m  m 


(12) 


for  values  of  k  between  0  and  /I.  The  range  in  over  which  aG  remains  monotonic 

is  determined  by  the  interaction  between  the  resonances  due  to  the  bulk  polar i tons 

and  the  tendency  of  the  polydispersity  to  smooth  out  such  resonance  effects.  Also, 

since  a  nonzero  imaginary  part  of  the  refractive  index  acts  to  damp  out  the 

resonances,  the  extinction  curves  for  «o(0)  display  the  most  abundant  resonant 

structure.  For  small  values  of  X  ,  (k )  increases  monotonically  with  respect  to  k; 

this  monotonic  behavior  of  a  (kf  continues  as  X  Increases  until  the  a  (0)  curve 

Dm  D 

crosses  the  nearest  neighboring  aD(k)  curve.  The  value  of  Xm  at  this  crossing  point 

defines  an  upper  limit,  Xmr  below  which  <*p(k)  is  monotonically  increasing  with 
respect  to  k.  The  value  of  X^  depends  both  on  the  polydispersity  of  the  aerosol 
and  on  the  value  of  the  real  part  of  the  complex  refractive  index.  The  table  shows 
the  approximate  values  of  X^  for  the  curves  of  figure  A-2.  Once  X^  becomes  larger 
than  Xmi’  the  behavior  of  a^(k)  is  not  generally  predictable  since  it  depends  in 
large  measure  on  the  resonant  structure  of  the  a_(0)  curve.  For  X  larger  than  X  . , 
it  does  not  appear  possible  to  determine  k  uniquely  from  extinction  data*  Figure 
A-3  shows  explicitly  the  variation  of  the  dimensionless  extinction  coefficient  with 
k  as  X^  increases  for  the  case  in  which  n  ■  2  and  o  -1.4.  The  curves  for  X  = 
1.38  and  X  -  1.44  bracket  X  .  for  this  case  (see  the^able). 
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Table.  Approximate  Values  of  the  Upper  Monotonicity  Limit 
for  the  Maaa  Median  Slae  Parameter 


1  % 

_ n _ _ _ I 

.1.33 

2.0 

3.0 

1.1 

3.55 

1.35 

0.89 

1.4 

3.55 

1.41 

0.81 

2.0 

4.68 

1.78 

0.89 

3.0 

9.77 

3.39 

1.23 

From  the  discussion  above.  It  may  be  concluded  that  the  extinction 
coefficient  will  be  monotonic  with  respect  to  k  for  the  range  of  optical  constants 
and  particle  sizes  likely  to  be  encountered  in  many  practical  problems  in  aerosol 
physics.  As  shown  in  figure  A-4,  it  now  becomes,  in  principle,  a  simple  to 
determine  k  numerically,  Is  the  extinction  coefficient  computed  from  equation  10 
using  the  current  estimate  for  k  and  the  experimentally  determined  size 
distribution;  a  is  the  measured  value  of  the  extinction  coefficient.  Successive 
estimates  of  k  are  made  until  |ac  -  ag|  <  6,  where  6  is  determined  by  the  precision 
to  which  the  extinction  coefficient  can  be  measured. 


2.3  Algorithm  for  computing  optical  constants. 

The  procedures  developed  above  may  be  combined  into  an  algorithm  for 
computing  the  optical  constants  as  shown  in  the  flow  chart  of  figure  A-5.  The 
computation  is  started  by  taking  n(A)  equal  to  a  constant  which  we  usually  choose  to 
be  1.3.  Successive  calculations  of  k(A)  and  n(A)  are  then  made  until  convergence  is 
obtained.  The  algorithm  was  implemented  by  means  of  a  FORTRAN  program  written  for 
use  on  a  Univac  1108  computer.  All  Mie  calculations  were  carried  out  by  means  of  a 
modified  version  of  Dave's  DBMIE  subroutine,4  and  the  original  version  of  the  SKK 
routine  used  in  this  study  was  developed  by  Querry.* 


3.  EXPERIMENTAL  PROCEDURE  AND  RESULTS 

An  extinction  spectrum  and  particle  size  distribution  for  o-phosphoric 
acid  were  measured  in  order  to  confirm  the  theory  developed  in  section  2.  Figure 
A-6  shows  a  diagram  of  the  experimental  arrangement  employed  for  these 
measurements.  The  o-phosphoric  acid  was  disseminated  in  a  22-m^  test  chamber  by 
spraying  the  acid  solution  through  a  pneumatic  nozzle;  the  acid  aerosol  was  stirred 
continuously  throughout  the  experiment  in  order  to  maintain  a  uniform  aerosol 
concentration.  An  Exotech  model  10-24  radiometer  with  a  circular  variable  filter 
monochrometer  was  used  to  scan  the  7-  to  ’4-pm  infrared  region  at  a  rate  of  15  scans 
per  minute;  the  path  length  through  the  aerosol  was  3.05  m.  The  radiometer  data 
were  recorded  on  analog  tape  with  a  Hewlett-Packard  model  3960  recorder.  Particle 


*  Querry,  M.  R.  Private  communication.  1967. 
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size  samples  were  taken  with  an  Andersen  model  2000  cascade  impactor,  and  the 
aerosol  mass  concentration  was  determined  from  samples  taken  on  glass  fiber  filters. 

A  Quad  Systems  model  721  digitizer  was  used  to  convert  the  analog 
radiometer  data  to  digital  form,  and  the  digitized  data  were  written  on  a  magnetic 
tape  for  use  In  subsequent  computer  processing.  The  extinction  spectrum  wa9 
determined  from  the  digitized  data  by  Beer's  law 

a  -  In  T  (13) 

e  CL 

where 

3 

C  -  aerosol  mass  concentration  (gm/m  ) 

L  -  optical  path  length  (m) 

T  •  transmittance 

The  measured  acid  concentration  of  the  aerosol  droplets  was  approximately  65%  by 
weigt  and  the  log-normal  particle  size  distribution  as  estimated  from  the  cascade 
impatLor  results  had  the  parameters  *  3.3  Um  and  a  *=  2.0.  Figure  A-7  compares 
the  experimentally  determined  extinction  spectrum  with  a  spectrum  computed  using  the 
optical  constant  data^  for  65%  by  weight  of  o-phosphoric  acid.  The  experimental 
spectrum  was  produced  by  averaging  45  radiometer  scans;  and  figure  A-8  shows  a  plot 
of  the  percent  difference  between  the  computed  and  experimental  spectra  as  a 
function  of  wavelength. 

4.  OPTICAL  CONSTANT  RESULTS 

The  experimental  data  for  o-phosphoric  acid  were  analyzed  according  to 
the  procedure  described  in  figure  A-5;  the  known  value  of  n(X)  was  taken  to  be  1.604 
at  10. 0  um.23  After  four  iterations,  the  calculations  converged  to  produce  the 
solution  shown  in  figure  A-9  for  n(X)  and  in  figure  A-10  for  k(X).  For  purposes  of 
comparison,  n(X)  and  k(X)  for  65%  o-phosphoric  acid  are  also  plotted  in  figures  A-9 
and  A-10;  figures  A-ll  and  A-12  are  plots  of  the  percent  difference  between  the 
computed  values  of  the  optical  constants  and  the  measured  values  for  65%  o- 
phosphoric  acid.  Notice  that  the  spectral  structure  of  the  percent  difference  for 
a(X)  (figure  A-8)  is  very  similar  to  the  spectral  structure  of  the  percent 
difference  for  k(X)  (figure  A-12).  This  suggests  that  differences  between  the 
computed  and  measured  values  of  the  optical  constants  are  due  to  experimental  error 
in  measuring  the  extinction  spectrum  and  particle  size  distribution.  In  order  to 
determine  the  amount  of  error  introduced  by  the  computational  procedure,  the 
computed  spectrum  from  figure  A-7  was  used  as  the  input  extinction  spectrum  from 
which  the  optical  constants  were  to  be  computed.  The  computation  converged  after 
five  iterations,  and  the  results  for  the  optical  constants  are  compared  with  the 
measured  65%  o-phosphoric  acid  data  in  figures  A-13  and  A-14.  The  average  percent 
difference  between  the  computed  and  the  measured  results  for  (n,k)  was  less  than  1 
percent  over  the  7-  to  14-pm  spectral  region.  Therefore,  it  was  concluded  that  the 
dominant  contribution  to  the  difference  between  the  optical  constants  computed  from 


the  o-phoephoric  add  extinction  data  and  the  neaaured  values  vas  due  to  error  in 
measuring  tha  extinction  spectrum  and  the  particle  else  distribution  of  the  aerosol. 

5.  CONCLUSION 

We  have  described  an  entirely  new  approach  to  the  determination  of  the 
optical  constants  of  particulate  materials.  This  new  method  avoids  the  experimental 
difficulties  and  errors  associated  with  the  traditional  methods  of  finding  the 
optical  constants  of  aerosol  materials.  It  is  an  iterative  procedure  which,  after 
an  initial  estimate  of  n(X)  over  the  wavelength  spectrum,  employs  a  simple  numerical 
algorithm  to  fit  the  k(X)  spectrum  to  a  measured  extinction  spectrum  and  the  KK 
dispersion  to  establish  a  new  estimate  of  the  n(X)  spectrum.  The  values  of  k(X)  are 
restricted  to  the  range  0  to  /?,  and  the  numerical  algorithm  for  finding  k(X)  is 
valid  for  a  wide  range  of  n(X)  and  particle  size  distributions  as  shown  in  the 
table.  The  value  of  n(X)  must  be  known  for  one  wavelength  in  the  spectral  region 
being  Investigated. 

The  theory  of  the  method  was  developed  In  detail  and  was  successfully 
applied  to  the  determination  of  the  optical  constants  of  an  o-phosphorlc  acid 
aerosol  in  the  7-  to  14-iim  Infrared.  The  numerical  procedure  was  shown  to  Introduce 
an  error  of  less  than  1  percent  in  the  determination  of  the  o-phosphoric  acid 
optical  constants. 
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Figure  A-l.  leopleths  in  n,k  Space  of  the  Dimensionless  Extinction 
Coefficient  for  Spheres  with  the  Indicated  Size  Parameters 


OiMIfiSiOfilfSS  txnwcno*  COiFFJCttM  OlMl  NSlONU  SS  IXtllWillON  CO»»»tci|*T  DIME  NSlOM  I SS  I  XT  INC  I  ION  COI  F  f  ICif  NT  OMSi«*SKM*lt«  IMTMtCTKM  COCFFKttttT 


Figure  A-2.  Dimensionless  Extinction  Coefficient,  ap,  for  Log-Normally 
Distributed  Spheres  as  a  Function  of  Mass  Median  Size  Parameter,  X 

IQ 

Individual  curves  are  for  k  ■  0,  0.2,  0.4,  0.6,  0.8,  1.0,  2  and  are 
easily  identified,  since  is  monotonically  increasing  with  respect  to  k  at  the 

left-hand  side  of  the  plots. 
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12.5 
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Figure  A-3.  Variation  of  the  Dimensionless  Extinction  Coefficient  with  k  as  Increases 
The  curves  for  Xm  -  1.38  and  Xm  “  1.44  bracket  X_o  • 
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Figure  A-4.  Illustration  of  the  Numerical  Procedure  for  Calculating 

a  Is  the  measured  extinction  coefficient  and  a  Is  the 
e  c 

extinction  coefficient. 


k 

computed 
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Index. 


Figure  A-5.  Flowchart  of  the  Algorithm  for  Determining  the 
Optical  Constanta  n,k  from  Extinction  Measurements 

nc  Is  an  Initial,  assumed  value  of  the  real  part  of  the  refractive 
nQ  Is  typically  taken  to  be  1.3. 
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flNFRARED 

SOURCE 


WR  ■  Radiometer  Window  -  Polyethylene 
Ws  -  Source  Window  -  Polyethylene 
L  -  Path  Length,  3.05  m 

a 

V  -  Chamber  Volume,  22  m 

Figure  A-6.  Experimental  Arrangement  Used  for  o-Phosphoric 
Acid  Aerosol  Measurement 
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Comparison  of  the  Experimentally  Determined  Extinction  Spectrum  (Dashed  Line) 
a  Computed  Extinction  Spectrum  (Solid  Line)  for  65%  o-Phosphorlc  Acid 


I 


Figure  A-8.  Percent  Difference  between  the  Computed  and 
Measured  Extinction  Spectrum  as  a  Function  of  Wavelength 


Comparison  of  n(X)  Determined  from  Extinction  Measurements  (Solid  Curve) 
h  the  Measured  Values  for  652  o-Phosphoric  Acid  (Dashed  Curve) 


Figure  A-10.  Comparison  of  k(A)  Determined  from  Extinction  Measurements  (Solid  Curve) 
with  the  Measured  Values  for  65%  o-Phosphoric  Acid  (Dashed  Curve) 


Perce**  Difference  between  n(X)  Determined  from  Extinction  Measurements 
the  Values  of  n(X)  for  65%  o-Phosphoric  Acid 


Con*>arison  of  n(X)  Derived  from  a  Spectrum  (Circles)  Computed  from  Lorenz-Mie  Theory 
with  the  Real  Parts  Used  in  the  Computation  (Solid  Curve) 


Comparison  of  k(A)  Derived  from  a  Spectrum  (Circles)  Computed  from  Lorenz-Mie  Theory 
with  the  Imaginary  Parts  Used  in  the  Computation  (Solid  Curve) 
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Figure  A-4.  Illustration  or  cne  Numerical  rroceouic  lul  ^ 

®e  i*  the  measured  extinction  coefficient  and  a 
extinction  coefficient.  c 
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